:arXiv^gr^qc/95K)7ibOvi" 29 Jul 1999: 



Semi-analytical Solution of Dirac equation in Schwarzschild Geometry 



Banibrata Mukhopadhyay 

and 

Sandip K. Chakrabarti 
Theoretical Astrophysics Group 
S. N. Bose National Centre For Basic Sciences, 
JD Block, Salt Lake, Sector-Ill, Calcutta-700091 

To appear in Classical and Quantum Gravity 

Separation of the Dirac equation in tiio spacotime around a Kerr black liole into radial and angular 
coordinates was done by Chandrasekhar in 1976. In the present paper, we solve the radial equations 
in a Schwarzschild geometry semi-analytically using Wentzel-Kramers-Brillouin approximation (in 
short WKB) method. Among other things, wo present analytical expression of the instantaneous 
reflection and transmission coefficients and the radial wave functions of the Dirac particles. Complete 
physical parameter space was divided into two parts depending on the height of the potential well 
and energy of the incoming waves. We show the general solution for these two regions. We also 
solve the equations by a Quantum Mechanical approach, in which the potential is approximated by 
a scries of steps and found that those two solutions agree. Wc compare solutions of different initial 
parameters and show how the properties of the scattered wave depend on these parameters. 
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I. INTRODUCTION 

The spacetime around an isolated black hole is flat and Minkowskian at a large distance where usual quantum 
mechanics is applicable, while the spacetime closer to the singularity is curved and no satisfactory quantum field 
theory could be developed as yet. However, occasionally, it is useful to look into an intermediate situation when a 
weak perturbation (due to, say, gravitational, electromagnetic or Dirac waves) originating from infinity scatters from 
a black hole. The resulting wave is partially transmitted into the black hole through the horizon and partially scatters 
oS^ from it to infinity. In the linearized ('test field') approximation this problem has been attacked in the past by 
several authors [1-4] . These methods are mostly numerical and most of the solutions obtained so far is for particles 
of integral spin only. Chandrasekhar [3-4] separated the Dirac equation in Kerr geometry into radial and angular 
parts. These works were extended to other spacetimes, such as in Kerr-Newman geometry [5], and around dyon black 
holes [6]. Subsequently, Chakrabarti [7] solved the angular part of the Dirac equation in Kerr geometry and gave 
the eigenvalues of the equation. These and the present works mostly concern scattering off tiny black holes and thus 
changing the incoming solution appreciably into an outgoing solution. Scattering effects from larger black holes could 
be studied by phase shift analysis and these has also been done recently [8]. 

In the present paper, we attack a simpler problem to have a 'feel' for the complete solution when the black hole 
is non-rotating. In the next Section, we present the basic equations. In §3, we classify the parameter space in terms 
of the physical and unphysical regions and present the method we adopt to solve the equations. In §4, we present 
a complete solution. In §5, we present solutions using a classical method in which the potential is approximated 
by a series of steps and then compare solutions of these two methods. In §6, we also compare solutions of various 
parameters and show how a Schwarzschild black hole distinguishes incoming particles of various masses. Finally, in 
§7, we draw our conclusions. 
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II. BASIC EQUATIONS OF THE PROBLEM 



Following Chandrasekhar [4] , the radial part of the Dirac equation is easily reduced into a Schrodinger like equation. 

However, because the spin-spin coiipling term is absent in the Schwarzschild geometry, the radial equation is much 
simpler to deal with. The eigenvalue of the angular equation for spin ±1/2 is trivially obtained as {1 + 1/2)'^ [7, 
9-10] where I is the orbital quantum number. In what follows, we choose 1 = 1/2 throughout for concreteness. This 
eigenvalue turns out to be the separation constant A of the original Dirac equation [4] . Here we solve the equation for 
one possible value of separation constant A (for / = 1/2, A is unity). In future we plan to explore the nature of the 
solutions for other orbital quantum numbers. 

Presently, we need to solve only the following coupled radial equations [4]: 

A5I?oi?_i = (1 +impr)A5i?_^i, (1) 
A^©tA^i?+i = (1 - impr)R_.. (2) 



where, 



^ ^ iK ^ (r-M) 
V^ = dr + -^+2n '- ^ 

A = - 2Mr, 



K = r'' 



u. 

Here n = integer, a = frequency of incoming Dirac wave, M = mass of the black hole, m = azimuthal quantum 
number, m,p = rest mass of the Dirac particle; p indicates particle, = radial wave function for spin up 

(down) particles. Vq is the complex conjugate operator. It is to be noted that the dimensionless unit is chosen, so 
that G = % = c = 1. The radial equation here is in coupled form. We can decouple it and express the equation 
either in terms of spin up or spin down wave function. However, it is more convenient to follow Chandrasekhar's [4] 
approach by which the basis was changed along with the independent variable r. That way, the coupled equation was 
reduced into two independent one dimensional wave equations since they are easier to solve. 
We first define 

r* = r + 2Mlog \r - 2M\ , (3) 

where, r > r+(= 2M), 

— = --, (4) 



and choose A^R,! = P,i, R_i = P_i. 

'2 ' 2 ' 2 2 

In terms of r*, the operators take the form: 



and 



We choose 9 = tan ^{wLpr) which yields. 



^\ , d 



cos^ = —77- 9~oT) sin^ = 



2 



and 



(1 ± irupr) = exp{±i^)^/{l + rripr'^). 
Following exactly Chandrasekhar's [4] approach we write 



and 



-P_ 1 = ip_iexp 



1 . , , 
— -itan [nipr) 



+-ztan [mpr) 



Finally, a choice of f* = r* + ^tan ^(mpr) yields df* = (-"^ + i+m'^r'^ ) ^'^*- 
With these definitions, the differential equations (1-2) are re-written as 



— -W]Z.= iaZ. 
an 



and 



where, Z± = tl>_^_i ±'tjj_i, and 



- + W]Z. = ^aZ,, 



W = 



A^(l + mlr'^f/'^ 



(5) 



(6) 



(7a) 



(76) 



(8) 



One important point to note: the transformation of spatial coordinate r to r* (and f*) is taken not only for math- 
ematical simplicity but also for a physical significance. When r is chosen as the radial coordinate, the decoupled 
equations for independent waves show diverging behaviour. However, by transforming those in terms of r* (and f*) 
we obtain well behaved functions. The horizon is shifted from r = r_|_ to = — oo. 

^Prom the above set of equations, we readily obtain a pair of independent one- dimensional wave-equations. 



+ a^]Z± = V±Z±, 



(9) 



where. 



dr* 



(10) 



A3(l + m2r2)3/2 
[r^{l + mlr^)+mpA/2a]- 



r[A^(l + m2r2)3/2 ± ((r - M)(l + m^r^) + Sw^rA)] 



AiCl + to2„2\5/2 

T .,.,^,,' ,3 [2r(l + my) + 2my + mp{r - M)/a]. 



(11) 



III. PARAMETER SPACE AND METHOD TO SOLVE EQUATIONS 



Wc obtain solutions by employing WKB [11-12] method and then imposing strict boundary conditions on the 
horizon, so that the reflection coefficient is zero and transmission coefficient is unity at the horizon. After establishing 
the general solution, we present here the solution of eq. (9) for two sets of parameters as illustrative examples. 
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It is advisable to choose the parameters in such a way that there is a significant interaction between the particle 
and the black hole. This is possible when the Compton wavelength of the incoming wave is of the same order as the 
Schwarzschild radius of the black hole, i.e., 



2GM n 



here, we are choosing G = ;i; = c= l, so 



Again, in the case of Schwarzschild geometry, frequency of the incoming particle (or, wave) will be of the same order 
as inverse of time. So, 



c3 



2GM 

Using the units as before, one can write, 

rrip ~ cr ~ (2M)"\ (12) 

In principle, however, one can choose any values of a and nip for a particular black hole and the corresponding solution 
is possible but we shall concentrate the region of parameter space where the solution is expected to be interesting as 
pointed out above, namely, region close to rUp = a. In Fig. la, we draw this line. The parameter space is spanned by 
the frequency a and the rest mass of the incoming particle m^. It is clear that 50% of total parameter space where 
a < nip is unphysical, and one need not study this region. Rest of the parameter space (cr > mp) is divided into two 
regions ~ I: E > Vm and II: E < Vm, where Vm is the maximum of the potential. While in Region I, the wave is 
locally sinusoidal because the wave number k is real for the entire range of f * . In Region II, on the other hand, the 
wave is decaying in some region when E < V, i.e., where the wave 'hits' the potential barrier and in the rest of the 
region, the wave is propagating. We shall show solutions in these two regions separately. In Rcgion-I whatever be 
the physical parameters, energy of the particle is always greater than the potential energy and WKB approximation 
is generally valid in the whole range (i.e. -^jf- « 1). In cases of Region-II, energy of the particle is always less 
than the maximum height of potential barrier. Thus, at two points (where, fc = 0) total energy matches with the 
potential energy and in the neighbour of those two points WKB approximate method is not valid. They have to be 
dealt separately. In Fig. lb, we show contours of constant Wmax =niax(p-^) for a given set {a,mp) of parameters. 
The labels show the actiial vahies of Wmax ■ Clearly, except for parameters very close to the boundary of Regions I 
and II, WKB approximation is safely valid for any value of f* . One has to employ different method (such as using 
Airy Functions, see below) to find solutions in this region. 

IV. THE COMPLETE SOLUTION 
Solutions of Region I 

In this region, for any set of parameter, energy of the particle is always greater than the corresponding potential 
energy. We first re-write equation (9) as, 

^ + (a^-y+)Z+=0. (13) 

This is nothing but Schrodinger equation corresponding to the total energy of the wave . This can be solved by 
regular WKB method [11-12]. Let, 

fc(n) = V(^'-V+), (14) 

u{f^)= I A:(f*)df, -|- constant. (15) 
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Here, k is the wavenumber of the incoming wave and u is the Eiconal. The solution of the equation (13) is, 

A+ A- 

= ——exp{iu) H — j-exp{—iu). (16) 

with 

Al+A'i= k. (17) 

In this case all along cr^ > V+ and also « k, so WKB approximation is generally valid in the whole region. 

The quantity -p-^ falls off rapidly with distance. Thus, WKB is strictly valid at long distance only. 

It is clear that a standard WKB solution where and A_ are kept constants throughout should not be accurate, 
since the physical inner boundary condition on the horizon must be that the reflected component is negligible there. 
Thus WKB approximation requires a slight modification in which the spatial dependence of A± is allowed. On the 
other hand, at a large distance, where WKB is strictly valid, and A- should tend to be constants, and hence their 
difference is also a constant: 

A+-A_= c. (18) 
Here, c is determined from the WKB solution at a large distance. This along with (17) gives. 



This spatial variation, strictly valid at large distances only, should not be extendible to the horizon without correcting 
for the inner boundary condition. These values are to be shifted by, say, A^-h respectively, so that on the horizon one 
obtains physical R and T. Wc first correct reflection coefficient on the horizon as follows: Let A-h be the value of 
A_ on the horizon (see, equation (19)), 

c , V[2fc(r+)-c2] 
^-' = -2+ 2 • 

It is appropriate to use A- — A- — A-h rather than A_ since A- vanishes at r = r+. 
Incorporating these conditions, the solution (16) becomes, 

A+ A+ 

Zj^ = —j-exp{iu) H -j-exp{—iu). (20) 

V9 VQ 



with the usual normalization condition 



Al + Al=q. (21) 



where, A+ = A+ — A+h- Here, q is to be determined self-consistently by equating the asymptotic behaviour of this 
reflection coefficient with that obtained using WKB method. This q in turn is used to compute A+ = A^ — A^h, and 
therefore the transmission coefflcient T from eq. (21). In this way, normalization of i? + T = 1 is assured. 

Normalization factor — > fc as f* ^ oo and the condition << q is found to be satisfied whenever << fc is 
satisfied. This is the essence of our modification of the WKB. In a true WKB, A± are constants and the normalization 
is with respect to a (almost) constant k. However, we are using it as if WKB is instantaneously valid everywhere. Our 
method may therefore be called 'Instantaneous' WKB approximation or IWKB for short. Using the new notations, 
the instantaneous values (i.e., local values) of the reflection and transmission coefficients are given by (see, eq. 20), 

A2 

i? = (22a) 



.42 

T=— t. (226) 
Q 

Determination of Aj^h is done by enforcing R obtained from eq. (22a) as the same as that obtained by actual WKB 
method at infinity. 
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To be concrete, we choose one set of parameters from Region I. (A large number of solutions are compared in §6 
below.) Here, total energy of the incoming particle is greater than the potential barrier height for all values of f*. 
We use mass of the black hole, M = 1; mass of the particle, rup = 0.8, orbital quantum number, I = ^, azimuthal 
quantum number, m = — i, frequency of the incoming wave, cr = 0.8. 

^From eq. (9) we observe that there are two wave equations for two potentials V+ and V- . The nature of potentials 
are shown in Fig. 2. It is clear that potentials V± are well behaved. They are monotonically decreasing as the 
particle approaches the black hole, and the total energy chosen in this case (cr^) is always higher compared to V±. 
For concreteness, we solve using potential V+. Similar procedure can be adopted using potential V- to compute Z_ 
and its form would be 

A' — A' A' — A' 

Z_ = + +^ exp{iu') - ~ . ^ -^ exp{-iu'). (20') 

Note the occurrence of the negative sign in front of the reflected wave. This is to satisfy the asymptotic property 
of the wave functions which must conserve the Wronskian [4] . Since the coefficients should not change sign between 
infinity and the horizon (as that would tantamount to having zero amplitude, i.e., unphysical, absence of either the 
forward or the backward component) the same sign convention is followed throughout the space. Local values of the 
reflection and transmission coefficients could also be calculated in the same manner. In the solution (eq. 20 and 20'), 
first term represents the incident wave and the second term represents the reflected wave. 

In Fig. 3 we show the nature of V+ (solid curve), k (dashed curve) and E{= a^) (short-dashed curve). The difference 
of E and V+ and therefore k goes up as the particle approaches the black hole. 

In Fig. 4, variation of 'local' reflection and transmission coefficients are shown. It is observed that as matter comes 
close to the black hole, the barrier height goes down. As a result, the penetration probability increases resulting in 
the rise of the transmission coefficients. At the same time, the reflection coefficient tends to be zero. It is to be noted, 
that, strictly speaking, the terms 'reflection' and 'transmission' coefficients are traditionally defined with respect to 
the asymptotic values. The spatial dependence that we show are to be interpreted as the instantaneous values. This 
is consistent with the spirit of IWKB approximation that we are using. 

The behaviour of the solutions with is not very different from what were shown in Figs. (3-4) except in a region 
very close to the black hole horizon where V+ and V- differs slightly (see. Fig. 2). 

Using the solutions of equations with potential V+ and V-, the radial wave functions R_^i and R_i for spin up 
and spin down particles respectively of the original Dirac equation are given below, 

-n (n a+cos{u — 6) -\- a-Cos{u + 6) a',cos{u' — 0) — a'_cos{u' + 6) 

Re I /iiA2 1 = — 1 — (23a) 

V ^ / 2\/fc 2^/k' 



Im 



. i\ a+sm(u — 9) — asiniu + 9) a',sm(u' — 9)+a'_sm(u' + 9) 



/ N a+cos{u + 9) + a^cos{u - 6) a'_^_cos{u' + 6) - a'_cos{u' - 6) 
Re I ii_ 1 1 = z= z= (2oc) 



Im 



/ ^ \ _ a+sin(u + 9) — a_sin(M — 9) a'_,_sin(M' + 9) + a'_sm{u' — 9) (2M) 



Here, a+ = (A+ — Aj^^)/ ^/{q/k) and a_ = (A_ — A^h) / \/\<l/k). Here, we have brought back k and k' so that these 
may resemble the original solution (eq. 16) using WKB approximation. and are the transmitted and reflected 
amplitudes respectively for the wave of corresponding potential K- . 

Figure 5 shows the resulting wave functions for both the spin +5 and spin — ^ particles respectively. The eiconals 
used in plotting these functions (see, crj. 23[a-d]) have been calculated by approximating V± in terms of polynomials 
(This was done since V± as presented in eq. 10 is not directly integrable.) and using the definition M(r*) = J \/{o''^ — 
V±)dr*. Note that the amplitude as well as the wavelength remain constants in regions where k is also a constant. As 
discussed before, the wave functions are almost sinusoidal close to the horizon and at a very large distance (albeit with 
different frequencies). Since the net current (|P+i |^ — |P_i P) is conserved, probability of spin +5 is complimentary 

to the probability of spin — ^ particles respectively. 
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Solutions of Region II 



Here we study the second region where for any set of physical parameter total energy of the incoming particle is 
less than the maximiim height of the potential barrier. Thus, the WKB approximation is not valid in the whole 
range of f*. In such regions, the solutions will be a linear combination of Airy functions because the potential is 
approximately linear in f* in those intervals. At the junctions one has to match the solutions with Airy functions 
along with the solution obtained by the WKB method. In the region where the WKB approximation is valid, local 
values of reflection and transmission coefficients and the wave functions can be calculated easily by following the same 
method described in Case I. In other regions, the equation reduces to 

^-XZ,=0, (24, 

where, a; = /Js (f* — p), /3 is chosen to be positive and p is the critical point where the total energy and potential 

energy are matching. 

Let Z^{x) = x^Y{x) and considering region .x > the equation (24) reduces to 

By making yet another transformation, 



^=lx"^, (26) 



we obtain. 



This is the modified Bessel equation. The solution of this equation is Z+i (^) and /_i (^). Hence, the solution of eq. 
(27) will be, 

Z+{x) = a;^[Ci7+i(0 + C2/_i(0]. (28) 
When a; < the corresponding equation is, 

which is the Bessel equation. The corresponding solution is 

Z+{x) = \x\^DiJ^i{0 + D2J-i{0], (30) 

where J± and I± are the Bessel functions and the modified Bessel functions of order i respectively. 
The Airy functions are defined as 

Ai{x) = ^x^I_.{0-I+^,m, ^>0, (31) 

Aiix)^^\x\i[J_^{0 + J+^m, x<0, (32) 

Bi{x) = ^xi[I_i{0 + I+UO], a;>0, (33) 



v/3 
1 

73' 



Bi{x) = ^\x\^J_,(0-J+M x<0. (34) 
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In terms of Airy functions, the solutions (28) and (30) can be written as 

Z+ = ^(C2 - Ci)Ai{x) + ^(C2 + Ci)Bi{x) for a; > 0, (35) 

Z+ ^1{D2 + Di)At{x) + ^{D2 - Di)Bi{x) for x < 0. (36) 

By matching boundary conditions it is easy to show that the solution corresponding a; > and that corresponding 
X < are continuous when Ci = —Di and C2 = -D2- 

To have an explicit solution, we choose the following set of parameters: M = 1, nip = 0.1, I = ^, m = and 
a = 0.15. 

In Fig. 6, we show the nature of V+ and V-. However, while solving, we use the equation containing V+ (eq. 9). 

Unlike the case in the previous Section, here cr^ is no longer greater than V± at all radii. As a result, fc^ may attain 
negative values in some region. In Fig. 7, nature of V+ (solid curve), parameter k (dashed curve) and energy E 
(short-dashed curve) are shown. Here, WKB approximation can be applied in regions other than f* ~ — 6 to — 1 and 
4 to 8 where k is close to zero and the condition << A; is not satisfied. In the region f * 8 to 4 around the 

turning point f* ~ 5.6088 the solutions turns out as [13] 

Z+ = 1.858386Ai(x) + 0.600610914Bi(a;). (37) 

Similarly, the solution from —1 to —6 i.e. around the turning point = —3.0675 can be calculated as [13] 

Z+ = 1.978145^i(a;) + 0.7168807Bi(a;). (38) 

It is to be noted that in the region ^ 4 to — 1, even though the potential energy dominates over the total energy, 
WKB approximation method is still valid. Here the solution will take the form "^"^^^"^ and ^^^7^^. Asymptotic values 

of the instantaneous reflection and the transmission coefficients (which are traditionally known as the 'reflection' and 
'transmission' coefficients respectively) are obtained from the WKB approximation. This yields the integral constant 
c as in Case I. /,From eq. 22(a-b) local reflection and transmission coefficients are calculated, behaviour of which are 
shown in Fig. 8. The constants A-h and A+h ai"© calculated as before. Note the decaying nature of the reflection 
coefficient inside the potential barrier. 



V. SOLUTION OF THE EQUATIONS BY STEP-POTENTIAL METHOD 

In the above sections we presented our semi-analytical solutions by WKB method with an appropriate boundary 
condition at the horizon. A numerical approach would be to replace the potential V^(r*) by a collection of step function 
as shown in Fig. 9a. Here, the solid steps approximate the dashed potential for rUp = 0.8 and a = 0.8. The standard 
junction conditions of the type, 

Z+,n = (39a) 

where, 

Z-\-^n ~ AT^SXp\ih'nT^^Trj\ -\- B^CXp[ ^Ari^*,n] 

and 



dZ+ dZ^ 



dr* df: 

where, 

dZ., 



\n+l 



(396) 



iknA>fiCXp(^iknV^^ri^ ikriB^expi^ iknT^^ri) 



at each of the n steps were used to connect solutions at successive steps. As before, we use the inner boundary 
condition, to be i? ^ at r* ^ —00. In reality, we used as many as 12000 steps to accurately follow the shape of the 
potential. Smaller step sizes were used whenever k varies faster. Fig. 9b shows the comparison of the instantaneous 
reflection coefficients in both the methods. The solid curve is from the WKB method of previous section and the 
dotted curve is from the step-potential method as we described here. The agreement is clearly excellent. 
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VI. BLACK HOLE: A MASS SPECTROGRAPH? 



In order to show that the black hole scatters incoming waves of different rest masses (mp) and of different energies 
(ct^) quite differently, wc show a collection of solutions in Figs. lO(a-d). In Figs. 10a, wc show reflection and 
transmission coeflScients for waves with parameters a = 0.8 (solid), 0.85 (dotted) and 0.90 (dashed) respectively with 
the same nip = 0.8. As the energy of the particle rises comparable to the height of the potential (which is solely 
dependent on nip at a large distance), the reflection coefficient goes down and the transmission coefficient goes up. 
In Fig. 10b, the real part of the wave Z^, corresponding to these three cases are shown. At f, = 0, the wave pattern 
is independent of a as the phase factor is trivially the same in all the cases. The dispersal of the wave with frequency 
is clear. Waves with smaller energy and longer wavelength are scattered with higher amplitude of Re(Z+) as the 
fraction of the refle cted wave goes up when energy is reduced. This behaviour is valid till R < 0.5 since the amplitude 
of Rc(Z+)= (1 + ^JTR)^/"^ . For R > 0.5, amplitude of Rc(2'+) goes down with energy. In Fig. lO(c-d), solutions are 
shown with varying the rest mass of the particles while keeping a fixed at 0.8. The solid curve, dotted curve and 
the dashed curves are for rup = 0.8, 0.76, 0.72 respectively. Most interesting aspect is that close to the black hole 
T-t ^ 0, the reflection and transmission coefficients as well as the nature of the wave are quite independent of the rest 
mass. This is understandable, as just outside the horizon, the potential is insensitive to nip. However, farther out, 
amplitude of Re(Z+) goes up as before when nip is raised as larger fraction of the wave is reflected from the potential 
(Fig.lOc). 

This interesting behaviour, for the first time shows that one could scatter a mixed wave off a black hole and each 
of the constituent wave would behave differently as in a prism or a mass spectrograph. 

VII. CONCLUDING REMARKS 

In this paper, we analytically studied scattering of spin-half particles from a Schwarzschild black hole. In particular, 
we presented the nature of the radial wave functions and the refiection and transmission coefficients. Our main 
motivation was to give an analytical expression of the solution which could be useful for further study of interactions 
among Hawking radiations, for instance. We verified that these analytical solutions were indeed correct by explicitly 
solving the same set of equations numerically using step-potential approach as described in Section V. Wc classified the 
entire parameter space in terms of the physical and unphysical regions and the physical region was further classified 
into two regions, depending on whether the particle 'hits' the potential barrier or not. We chose one illustrative 
example in each of the regions. We emphasize that the most 'interesting' region to study would be close to nip ^ a. 
However we pointed out (Fig. lb) that for nip < 0.3, WKB solutions cannot be trusted, and other methods (such as 
those using Airy functions) must be employed. 

We used the well known WKB approximation method as well as the step-potential method of quantum mechanics 
to obtain the spatial dependence of the coefiicients of the wave function. This in turn, allowed us to determine the 
refiection and transmission coefficients and the nature of wave functions. The usual WKB method with constant 
coefficients and (almost) constant wave number k is successfully applied even when the coefficients and wave number 
are not constant everywhere. Solution from this 'instantaneous' WKB (IWKB for short) method agrees fully with 
that obtained from a purely classical numerical method where the potential is replaced by a collection of steps (simple 
quantum mechanical approach). The resulting wave forms satisfy the inner and the outer boundary conditions. Our 
method of obtaining solutions should be valid for any black hole geometry which are asymptotically flat so that radial 
waves could be used at a large distance. This way wc ensure that the analytical solution is closer to the exact solution. 
In Region II, in some regions, WKB method cannot be applied and hence Airy function approach or our step-potential 
approach could be used. 

In the literature, reflection and transmission coefficients are defined at a single point. These definitions are mean- 
ingful only if the potential sharply changes in a small region while studies are made from a large distance from it. In 
the present case, however, the potential changes over a large distancx; and we are studying in these regions as well. 
Although we used the words 'reflection' and 'transmission' coefficients, in this paper very loosely, our definitions are 
very rigorous and well defined. These quantities are simply the instantaneous values. In our belief these are more 
physical. The problem at hand is very similar to the problem of reflection and transmission of acoustic waves from a 
strucked string of non-constant density where reflection and transmission occurs at each point. 

Unlike in the case of a Kerr black hole, the solution of the angular equation around a Schwarzschild black hole is 
independent of the azimuthal or meridional angles [5-7]. This is expected because of symmetry of the space-time. 
However, otherwise, the nature of the reflection and transmission coefficients were found to strongly distinguish 
solutions of different rest masses and different energies as illustrated in Figs. lO(a-d). For instance, when the energy 
of the wave is increased for a given mass of the particle, reflected component goes down. In regions where R > 0.5, 
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Rc(Z_i_) goes down with energy, but where R < 0.5, Re(Z-|_) goes up with energy. In any case, the wave scattered off 
are distinctly different. In a way, therefore, black holes can act as a mass spectrograph! For instance a mixture of 
waves should be splitted into its components by the black hole. Our method is quite general and should be used to 
study outgoing Hawking radiations also. This is beyond the scope of this paper and would be dealt with in future. 
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FIGURE CAPTIONS 



Fig. la: Classification of the parameter space in terms of the energy and rest mass of the particles. The physical 
region a > rUp is further classified in terms of whether the particle actually 'hits' the barrier or not. 

Fig. lb: Contours of constant Wmax= iiiax(-p-^) are shown to indicate that generally w « 1 and therefore WKB 
approximation is valid in most of the physical region. Labels indicate values of w. 

Fig. 2: Behaviour of potentials V+ (solid curve) and V- (dashed curve), as a function of f*. The parameters are 
chosen from Region I of Fig. 1. 

Fig. 3: Behaviour of V+ (solid curve), k (dashed curve), total energy E (short-dashed curve), as functions of r*. 

Fig. 4: Behaviour of local transmission (T, solid curve) and reflection {R, dashed curve) coefficients, as functions of 
f*. Close to the horizon, transmission is a hundred percent and reflection is almost zero. 

Fig. 5: Behaviour of (a) Re{Ry2A^/^), (b) Im(i?i/2Ai/2) , (c) Re(i?_i/2), (d) Im(i?_i/2), as a function of f*. 
Note the complimentary nature wave functions of the spin +i and spin — ^ particles. This is because the current is 
conserved. 
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Fig. 6: Behaviour of V+ (solid curve) and V- (dashed curve), as a function of r*. The parameters are chosen from 
Region II of Fig. l(a-b). 



Fig. 7: Descriptions are same as in Fig. 3. See text for the choice of parameters. 



Fig. 8: Descriptions are same as in Fig. 4. See text for the choice of parameters. 



Fig. 9a: Steps (soUd) approximating a potential (dotted) thus reducing the problem to that of a quantum mechanics. 
The parameters are rup = 0.8 and a = 0.8. 

Fig. 9b: Comparison of variation of instantaneous reflection coefScient R with the radial coordinate f* using analj^tical 
WKB method (solid) and numerical step-potential method (dotted). The parameters are rUp = 0.8 and a = 0.8. 

Fig. lO(a-d): Comparison of (a) reflection and transmission coefficients and (b) real amplitude of the wave function 
for nip ~ 0.8 and a ~ 0.80 (solid), 0.85 (dotted) and 0.90 (dashed) respectively, (c-d) Similar quantities for 
nip = 0.80, (solid) 0.76 (dotted) and 0.72 (dashed) respectively keeping a = 0.8 fixed. Higher reflection component 
enhances the wave amplitude, thus differentiating the incoming waves very clearly. 
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